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ABSTRACT 



This thesis presents a discussion of the problems 
involved in estimation of Circular Error Probable (CEP) . 
Several estimators are compared through simulation, under 
two models; the power model and the Rayleigh model. Several 
measures of effectiveness are calculated for each of the 
competing estimators. It is found that maximum likelihood 
estimation based on the power distribution performs well 
for "heavy tailed" distributions; the Rayleigh unbiased 
estimator performs well for most other situations. 
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INTRODUCTION 



I . 

Evaluation of a weapons system is frequently an expen- 
sive and time consuming task. Only relatively few tests 
may be allowed because of the high cost of testing. It is 
possible that a good weapon system could be rejected because 
of inefficient utilization of the small amount of data 
available. Hence, efficient use of the limited data is 
required to reduce this risk. 

Several authors have considered the measurement of 
weapon delivery accuracy using Circular Error Probable 
(CEP) . This measure of accuracy is defined to be the radius 
of a circle centered at the target which on the average con- 
tains fifty percent of impact points, that is, the median 
radial miss distance. A functional relationship between 
CEP and the distribution parameters is given by the 
following integral: 

CEP 

F (CEP) = / f (r) dr = 0.5 , 

0 

where f(r) is the probability density of radial miss 
distance . 

Various parametric models for the distribution of impact 
points about the target have been used in connection with 
such estimations. A discussion of these models is given by 
Jordan [1] . But, deficiencies in the models have been 
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noted. Specifically, they tend to provide poor fit to the 
"upper tail" of the radial distributions. That is, they 
provide less than satisfactory explanation of impacts ob- 
served far from the target. A commonly used model is the 
"Rayleigh model", in which impacts in the target plane are 
assumed to be distributed as the unbiased circular normal 
distribution. One might anticipate that Rayleigh-based 
CEP estimators might not perform adequately with a large 
proportion of actual test data. Several potential competitors 
to Rayleigh CEP estimators were therefore developed for 
comparison . 

Barr and Jordan [2] developed and proposed the power 
distribution model which is based on an infinite mixture of 
Rayleigh radial miss distance distributions. This distri- 
bution has properties compatible with the interpretation of 
o in the underlying Rayleigh distribution model; they sug- 
gested the use of maximum likelihood estimates of the param- 
eters of the power distribution which in turn provide a 
maximum likelihood estimate (MLE) of CEP. 

The jackknife is a technique which was proposed by 
Quenouille [10] for reducing bias in estimation problems. 

We tried the jackknife procedure to estimate CEP under the 
Rayleigh model. 

The sample median was also used to estimate CEP. Since 
CEP is defined as the median radial miss distance, a natural 
estimator for CEP is the sample median. 
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To measure the performance of the estimators, a variety 
of data sets from known distributions, such as uniform, 
Rayleigh, power and Weibull were generated by simulation. 
Samples from the known distributions were generated to 
test the capability of the estimators to handle widely 
disparate situations. 
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II. MAXIMUM LIKELIHOOD ESTIMATION BASED 
ON THE POWER DISTRIBUTION 

Barr and Jordan [1], [2] have given a discussion of 

CEP accuracy measurement under a power distribution model. 

They showed that estimating the parameters of the power 

distribution by the method of moments may not be adequate; 

moreover it is possible only when the ratio of the squared 

—2 2 

sample mean and variance r /s^ does not exceed tt/ 4 - if. 
Consequently, they suggested use of the maximum likelihood 
(M.L.) method. 

Analytically closed expressions for the M.L. estimators 
cannot be obtained, the solution of normal equations yielding 
the parameter estimates requires numerical techniques. In 
this section we examine a numerical method to obtain the 
M.L. estimates of the parameters of the power distribution, 
which in turn provide a M.L. estimate of CEP. These estimation 
procedures require use of a computer. 

A. MATHEMATICAL FORMULATION 

The probability density function of the radial miss 
distance R having a power distribution with parameters k and 
C is 



f R (r ; k, O 



2<r 

C 



(- 

C 



+ r 



K+l 



r > 0 , k , £ > 0 
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A C.D.F. for radial miss distance is 



F R (r ; k, O 



1 - (- 



+ 




< 



r>0, k, £ > 0 (1) 



The CEP is given by 

1 1 

CEP = F R _1 (|) = [ ( 2 K - l)^] 2 



Because radial miss distance is always non-negative, 
one can equivalently consider squared radial miss distance. 

In the case of the power model, this leads to some notational 

. . . 2 
simplification. Let Y denote the squared miss distance R . 

Then from (1) , 

F y ( y ) = 1 - (1 + Z)~ K ; y > 0 , K ,c > 0 

( 2 ) 

f y (y) = |(1 + ^)~ K ~ 1 ; y > 0 , k,^ > 0 

Here 5 and k are parameters to be estimated from observed 
impact data. The relationship between CEP and these 
parameters is 

1 

CEP 2 = F y _1 (|) = [ ( 2 K - 1)^] 



The method of maximum likelihood is based upon the 
likelihood function 4>(y) which is defined as the joint 
density function of a sample of size n; i.e.. 
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n n n y . _ 

<P (y) = tt f (y.) = —■ tt (i + -p)“ K “ 
i=l * 1 c i=l ^ 

Since it is easier to deal with suras rather than products 
when maximizing, c{> ( y ) is transformed into 

L (y) = log <p (y) 

n 

= n log k + nxlog £ - (k+ 1) 1 log(£+y.) 

i=l 1 

There is no loss of generality here because the maximum of 
a positive function occurs at the same point as the maximum 
of the logarithm of the function. The normal equations are 
then 



m(< , £) 



_ 9L(y) 

3 vc 



n 

- + nlog C ~ 2 logU+y .) 

K i=l 1 



0 



and 



g (<,0 



9L(y) 

U 



n< 

C 



(k+1) 



n 

E 

i=l 



1 



c + y. 



o 



The solution of m(£,£) =0 and g(K,£) = 0 gives the 

/\ 

parameter M.L. estimators k, t; . As mentioned before, in 
the present case, this requires numerical techniques. 



B. NUMERICAL SOLUTIONS: NEWTON'S METHOD 

This method is a Newton's iteration technique and the 
concept and procedure are discussed in [3], Newton's method 
provides a statistical criterion for stopping the iterations. 
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The pure form of Newton's method is given by the 



following iteration scheme: 



/\ _ I ^ m 

X j+1 = Xj - [F (X j ) ] X Vf (Xj) 



where 



_ A A _ 

j ’ ‘"j-'j 1 ' 



£ (x . ) = JililLj 



3k 



3; 



First we compute 



Let 



3L(y) : 

A 

a< 



n 



n 



x + n log 5 - Z log U+y ■ ) 

K i=l 1 



3L(y) = 
3 C 



n< 



n 



- (k+1) Z 



i=l c + Yi 



9 2 L(y) 

9k 



-n 

-2 

K 



3 2 L(y) = n _ £ 1 

A A A A 

3<3 C C i=l C + Yi 



3 2 L(y) 

-2 

3C 



- ~ n < 



n 



+ (k+1) Z 



i=i U + y L ) 



D = 3 2 L(y) . 3 2 L(y) _ ^ 3 2 L (y ) ^ 2 

^ ^ 2 ^ A 

3k 3; 3k3c 
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If D were identically equal to zero, this method would not 
be applicable. Now take 



-1 ~ 

F (X) = 



3 2 L (y) _ 9 2 L(y) 

A a 

9C 9k9£ 

3 2 L(y) 3 2 L (y ) 

A A AO 

9k9C 9k 



D 



In order to develop the termination criterion, consider 
the statistic for (for reference see [4]) 

B = Vf (X.)F _1 (X.) Vf (X.) _1 
J 1 3 

B is asymptotically a Chi-square r. v. with 2 degrees of 

freedom, therefore it might be reasonable to continue the 

iteration until the statistic B becomes less than the 100a 

percentile point of the chi-square distribution with two 

2 

degrees of freedom, x ( a ) • 
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III. MAXIMUM LIKELIHOOD ESTIMATION BASED 
ON THE RAYELIGH DISTRIBUTION MODEL 



Suppose impacts, measured in a suitable plane, have 

range (X) and deflection (Y) components of impact that are 

jointly uncorrected bivariate normal distributed with mean 

2 

at the target and common unknown variance a in both direc- 
tions. That is, for the present development, it will be 
assumed that the joint density of miss distances about the 
target is given by 

,, , 2,-1 - (x 2 +y 2 ) /2a 2 

f (x,y) = (2tto ) e 2 , -» < x,y < 

It is easy to derive the density of radial miss distance. 
1/2 

Let R = (X+Y) and transform to polar coordinates to 
obtain 

, R 2 

1 2tt r 2 

F (r) = P [ R <r ] = — / / Re 20 dR d0 

2tto 0 0 

- 2 

r 

9 2 o 

= 1 - e 0 ; r >_ 0, o > 0 . 

This distribution of the radial miss distance is referred 
to as either the radial normal distribution or the Rayleigh 
distribution . 

Then the density function of the radial miss distance, 
f R (r) is 
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f R<r> 



r 2g‘ 

~~7 e 

o 



r > 0 , a > 0 



To find the relationship between CEP and a, we calculate 



CEP = F r 1 (i) = J2 In 2 a : 1.1774 o , 



which is a well known relation. Hence, if we want to estimate 
CEP with this model, first we might estimate the population 
standard deviation. The maximum likelihood estimator for 
a is easily found to be 



a 




n 2 2 
Z r-V 

i=l 1 



/ 



where r^, ..., is the sample of observed radial miss 

distances. This estimator has a slight bias, however. An 
adjusted unbiased estimator for CEP is 



CEP = 1.177 4 (/n SlUl ) g, (3) 

r(|(2n+l) ) 

where the quantity in parenthesis is an unbiasing factor 

A 

for o. It has been shown by Chapman and Robbins [6] that 
CEP = 1.1848 a for sample size 20 is a best estimator of CEP 
under the normality and other assumptions stated above. 
However other estimators, while not as efficient as the best 
estimator, may have more robustness under variations in 
the assumed distribution of radial miss distance. 
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The term "jackknife" has been used by Miller [9] and 

Quenouille [10] to describe a method of modifying a biased 

estimator to bring it closer to an unbiased condition. 

From (3) we see the M.L. CEP estimator has a bias term of 

order — . We are thus led to try to examine the jackknife 
/n 

method to reduce this bias in the estimator. Brillinger [7 ] 
has provided a discussion of the asymptotic properties of 
estimators obtained with the jackknife method when it is 
applied to M.L.E's. 

The jackknife procedure is based on dividing the data 
into groups, obtaining estimates from combinations of the 
groups, and then averaging these estimates. We apply this 
to our problem, as follows. The maximum likelihood estimator 
for c is 



-2 

a 




r . 
D 



2 



Divide these data into m groups, where r^, ..., r^ is the 

sample of n observed radial miss distances and the size of 

~2 

each subgroup is k. First estimate the parameter o based 
on all data 



A 




2 



1 n 
-t=- lx.' 

2n j=l ^ 



~2 th 

next estimate o based on all data except those in the l 

group 
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-2 1 5 2 , 0 

°-i = 2H .5. r j ; 1 = X ' 2 ' • 



• t m. 



Then compute the pseudo values 



~ 2 ~ 2 ~ 2 
a. = me - (m-l)a . ; i = 1, 2, 

1 o -1 ' 



. , m. 



Finally, compute the average of the pseudo values, 

m ^ o 

Z a. 

e 2 = i =1 1 



m 



and estimate the CEP by 






CEP = i/2 In 2 a 
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IV. THE SAMPLE MEDIAN 



Since CEP is defined as the median radial miss distance, 
a natural estimator for CEP is the sample median (m) . By 
definition; the median is the 5C) percentile of the distri- 
bution, i . e . , 

P(R £ m) 21 0.50 and P(R >. m) >_ .50 

Since F (CEP ) = 0.50 and R is assumed to have a continuous 
distribution , 



~ -11 
CEP = m = F (i) 

Suppose that the n radial miss distance observations of 
the sample are arranged in order of magnitude, r^, r (2) ' 
..., r ( n )- We shall consider only the estimates of the 
population median defined as follows: 

(a) If the sample size is odd, and n = 2p+l, the median 

t h 

(m) is taken as the (p+1) c value. 

(b) If the sample size is even, and n = 2p, the median 

i.r 

(m) is taken as the midpoint between the p 1 and 
(p+l)^ values or 



1 

7 



(r 



(P) 



+ 



r (P+1) 



) . 
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The distribution of the median will tend to be symmetri- 
cal, and should be very nearly normal in almost all popula- 
tions. The distribution of the median in the case of small 
samples from a normal population, whether the sample size 
is odd or even, will tend to be normal [Hoji [11]) . The 
normal form is approached very rapidly. The value of the 
standard deviations of the median for n even and n odd 
approach the same limit. 



where a is the standard deviation in the population sampled. 

We use the order statistics to examine the sample median 
under the Rayleigh distribution model (see Inselmann [12]), 
and find its bias and variance so as to enable one to make 
a comparison with other estimates. If n is even, the dis- 
tribution of sample median is given by the following: 




g ( r ( o. 5 n ) ) = { n ! / (n- [ 0 . 5n] ) ! ( [0.5n]-l) !}f(r)F(r [0>5n] 



[0. 5n]-l 




j j n- [ 0 . 5n] 



For convenience we let a = 0.5n, and 



C = n!/ (n- [0 . 5n] ) ! ( [0. 5n] -1) ! = (2a) './a \ (a-1) ! 
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Then we may write the distri- 



so the sample median is r . 
bution of sample median under the Rayleigh distribution model, 



2 2 , 2 
g(ra) = c exp ] [1-exp ] a ± [exp( — 

o 2 a 2 a 2a 

2 2 
r , , r (a+l)., r , , r . .a-1 

= c—2- [exp (- - — 2 — )H1 - exp ( 

o 2a 2a 



We will now consider computation of the mean of ra. 



oo 2 2, ,, 2 _ , 

E (ra) = c / ^ r [exp(- r ta ■■)] [ 1 -exp (- — s-)] 3 dr 

0 a 2a 2a 



” r 2 r , r 2 (a+lK, a ^ , ..K.a- 1 . r , r 2 ,,K, 

= c / — [exp (- ^ )1 Z (-1) ( K ) [exp( ^) ] dr 

0 a 2 a K =0 2 o 



a „ 1 , . . K , a- 1 . “ r 2 , r 2 (a+l+K) . . 

= c Z (-1) ( „ ) / -j exp ( 5 -) dr 

K =0 0 a 2 o 



and make the following transformation 
y = r 2 (a+l+K) / 2 a 2 , 
r = a J 2 y/ (a+l+K) , 

dy = [r (a+l+K) /a 2 ] dr . 
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To obtain, 



E(ra) = ca/T Z (-1) K f 3 " 1 ) / (a+l+K) 2 / J~y e y dy 

K=0 * 0 



3 

2 



,3, ^ . a; 1 i-dVk 1 ' 



= cor(|) -/2 * i 



K=0 (a+l+K) 3/2 



We can adjust the bias term and using the same transformation 
as above we may now find the second moment 



E< r a 2 ) 



a-1 



K ,a-l, 



K 



= c 



/ t^xp(- r ■ ( r 1) ) £ < ' 1) ( K > texpl--^)] dr 



0 0 



2o 



K=0 



2o 



2 a-1 (-l)V" 1 ) 

= 2C0 E sr 



K=0 (a+l+K) 

The variances of this sample median can now be computed 
using the well-known formula 



Var (r a ) = E(r a 2 ) " t E (r a ) ]2 * 
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V. SIMULATION METHODOLOGY 



The previous sections have been concerned with a 
development of different estimators for the median of the 
distribution of radial miss distances obtained under a set 
of fixed impact conditions. An analytical comparison is 
difficult to perform for these estimators, so computer 
simulation was used. This section summarizes the different 
models and estimators discussed in the previous sections, 
and includes an analysis of the results obtained from simu- 
lations. Although the sample problems do not represent actual 
weapon test results, an attempt has been made to generate 
data from a wide variety of "realistic" distributions. 
Therefore this simulation analysis should show relationships 
between our estimators of CEP, and indicate their relative 
robustness, and hence potential for use in a variety of 
applications. The alternative estimators are as follows: 

A 

CEP 1: maximum likelihood unbiased Rayleigh estimator 

A 

CEP 2: maximum likelihood Jackknife Rayleigh estimator, 

with subgroup sample size k = 1 

A 

CEP 3: maximum likelihood power estimator under Newton's 

method 

A 

CEP 4: sample median estimator 

For the purpose of our simulation comparison of these 
four CEP estimators, we generated samples of radial miss 
distances (or, in some cases, squared radial miss distance) 
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from several parent distributions: Rayleigh, uniform, power 

and Weibull distributions. Throughout, we assumed a sample 
size of 20 impact observations. The method of generation 
was by C.D.F. inversion of generated uniform (0,1) values. 
For each parent distribution (hereafter called a "case") 
we replicated estimation of CEP's 500 times. Summaries of 
the sample distributions of 500 values of 5 CEP's were 
developed . 

Eight cases were considered; in each, parameter values 
were used which gave a true CEP of 100. The cases are as 
follows (U is a uniform (0,1) random variable throughout): 
Case I. (Rayleigh distribution) 

2 

F ^ (y) = 1 - e y//2 ° with o 2 = 10-^/21n2 

CEP = J~2 In 2 a 

Generated value: Y = -2a^ In U = -14426.95 In U 

Case II. (Uniform distribution) 

f r (y ) = 2^0 y ' with y e (0 ' 200) • 

CEP = 100 

Generated value Y = 200 U . 
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Case III. (Power distribution) 



F R 2 (y) = 1 - (1 + |)" K 



y > 0 



> 0 



with k = 1 and £ = 10 ^ 

1 1 
CEP = U (2 K - 1) ] 2 



_1 

Generated value Y = (U K - 1) z, = 10^ - 1) 



Case IV. (Translated Rayleigh distribution) 
F 2 (y) =■ 1 - e " (Y-100) /2o 2 

K 

where 2a 2 = 14282.6809 

1 

CEP = [100 + 2a 2 In 2] 2 



2 

Generated value Y = - 2o In U + 100 



Cases V-VIII. (Weibull distributions) 

w 

v 

2 _y 
F R ( y ) = 1 - e 

_1_ 

CEP = (v In 2 ) 2w 



Generated value Y in the table below 
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I 






Four sets of parameters (w,v) were selected to give 
CEP = 100 and to provide distributions with varying upper 
tail thickness. One case (with w = 1) gives an exponential 
distribution with mean v, and this corresponds to the Rayleigh 
case. Case I above may be considered also as a Weibull case. 
The parameter pairs used, in order of increasing upper tail 
size, are given in the following table. 



1 



Case 


w 


V 


Y = 


(-vlnU) V ' 


V 


2 


10 8 /ln2 


Y = 


(-144260504.1 In U) 

2 


VI 


1.5 


10 6 /ln2 


Y = 


(-1442695.04 In U) 3 


I 


1.0 


10 4 /ln2 


Y = 


-4426.95 In U 


VII 


.5 


10 2 /ln2 


Y = 


(-144.27 In U) 2 


VIII 


.1 


10 * 4 /ln2 


Y = 


(-3.623886 In U) 10 
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VI. SIMULATION RESULTS 



As mentioned above, we tested four CEP estimators from 
eight different cases. Major comparison measures include 
bias, variance and mean square error. Others are included 
in the tables that follow. 
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CASE I 



RAYLEIGH DISTRIBUTION 



Y = - 14426.95 In U 



* 



MOE 


CEP 1* 


CEP 2 


CEP 3 


CEP 4 


MEAN 


99.93 


99.94 


95.96 


107 


STD DEV 


11.66 


11.68 


12.22 


38.8 


MSE 


136.03 


136.37 


165.77 


1549 


GED MEAN 


99.25 


99.25 


95.17 


100 


VARIANCE 


136.03 


137.37 


149.45 


1500 


COEFF VAR 


0.12 


0.12 


0.13 


0.36 


SKEWNESS 


0.15 


0.14 


0.049 


0.35 


KURTOSIS 


0.10 


0.085 


0.307 


0.34 


RANGE 


71.38 


70.87 


83 ..60 


252 


MINIMUM 


67.25 


67.29 


53.74 


17.7 


.10 QUANTILE 


84 . 82 


84.80 


80.99 


58.4 


.25 


92.12 


92.09 


87.67 


82.5 


.50 


99.63 


99.63 


95.96 


105 


.75 


107.49 


107.34 


104.00 


132 


.90 


114.83 


114.83 


111.32 


156 


MAXIMUM 


138.63 


138.16 


137.34 


269 



CEP 1: maximum likelihood unbiased Rayleigh estimator. 

A 

CEP 2: maximum likelihood Jackknife Rayleigh estimator, 

with subgroup sample size k = 1. 

/N 

CEP 3: maximum likelihood power estimator under 

Newton's Second method. 



CEP 4: sample median estimator. 
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CASE II 



UNIFORM DISTRIBUTION 





Y 

/V 


= (200 U) 

/V 


A 


A 


MOE 


CEP 1 


CEP 2 


CEP 3 


CEP 4 


MEAN 


98.61 


96.52 


91.24 


102 


STD DEV 


9.87 


9.73 


17.08 


41.7 


MSE 


109.32 


107.03 


368.63 


1744 . 0 


GED MEAN 


96.09 


96.01 


* 


91.5 


VARIANCE 


97.83 


94.92 


291.89 


1740 


COEFF VAR 


0.10 


0.10 


0.19 


0.40 


SKEWNESS 


- 0.21 


- 0.20 


- 1.83 


0.005 


KURTOSIS 


0.61 


0.60 


5.54 


- 0.56 


RANGE 


67.28 


66.35 


125.60 


191 


MINIMUM 


59.52 


59.97 


0.0 


3.75 


.10 QUANTILE 


84.19 


84.25 


72.1 


45.7 


. 25 


90.10 


90.13 


85.30 


742 


.50 


96.84 


96.74 


94.26 


102 


.75 


102.96 


102.75 


101.52 


132 


.90 


108.66 


108.40 


107.62 


158 


MAXIMUM 


126.80 


126.32 


125.60 


194 


/s 

* some CEP = 


0.0 
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CASE III 



POWER DISTRIBUTION 



Y = 10 4 • (i - 1) 



MOE 


CEP 1 


CEP 2 


CEP 3 


CEP 4 


MEAN 


209.61 


229.33 


100.71 


113 


STD DEV 


190.15 


257.90 


21.07 


46.0 


MSE 


2460 


83000 


444.36 


2279 


GED MEAN 


177.76 


185.62 


98.60 


10 3 


VARIANCE 


3615.71 


66515.7 


443.86 


2110 


COEF VAR 


0.90 


1.12 


0.21 


0.40 


SKEWNESS 


6.39 


7.32 


0.84 


0.76 


KURTOSIS 


57.54 


72.76 


2.05 


1.48 


RANGE 


2264.4 


3220.29 


158.52 


337 


MINIMUM 


74.59 


75.08 


4 8.32 


17.2 


.10 QUANTILE 


108.29 


108.79 


76.88 


57.1 


. 25 


126.07 


127.96 


86.19 


82.5 


. 50 


161.63 


166.20 


99.25 


107 


,75 


220.09 


228.03 


112.71 


140 


.90 


331.77 


364.51 


124.77 


172 


MAXIMUM 


2339.01 


3295.38 


206.84 


354 
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CASE IV. TRANSLATED RAYLEIGH DISTRIBUTION 



Y = - 14282. 



MOE 


CEP 1 


MEAN 


99.79 


STD DEV 


11.56 


MSE 


133.74 


GED MEAN 


99.12 


VARIANCE 


133.70 


COEF VAR 


0.12 


SKEWNESS 


0.15 


KURTOSIS 


0.10 


RANGE 


70.75 


MINIMUM 


67.44 


.10 QUANTILE 


84.81 


.25 


92.04 


.50 


99.48 


. 75 


107.28 


.90 


114.56 


MAXIMUM 


138.19 



In U + 100 



A A 



2 


CEP 3 


CEP 4 


79 


95.85 


108 


58 


12.64 


38.3 


09 


176.96 


1524 


11 


94 . 78 


100 


05 


159.74 


1460 


12 


0.13 


0.35 


15 


- 0.56 


0.36 


08 


4.59 


0.34 


25 


135.40 


248 


46 


7.73 


20.4 


78 


81.18 


59.0 


00 


87.60 


82.7 


48 


96.18 


105 


12 


103.97 


132 


59 


110.99 


156 


71 


143.13 


268 



6809 

A 

CEP 

99 

11 

134 

99 

134 

0 

0 

0 

70 

67 

84 

92 

99 

107 

114 

137 
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CASE V. WEI HULL DISTRIBUTION w 2 

I 

Y » (- 1442 0050*1 . I In U) ~ 



)E 


CEP 1 


CEP 2 


CEP 3 


CEP 


:an 


86.18 


85.80 


S '! .25 


100 


'D DEV 


5.23 


5.19 


11.47 


19. 


;e 


218.34 


22S.61 


379.57 


382. 


:d mean 


86.03 


8 5.64 




98. 


\RIANCE 


27.35 


26.97 


131.51 


382 


OEF VAR 


0.06 


0.06 


0.14 


0 . 


KEWNESS 


0.029 


0.028 


- 5.47 


- 0 . 


1 IRTOSIS 


0.28 


0 . 27 


38.40 


0 . 


vNGE 


32.67 


32.46 


112.79 


121 


NIMUM 


71.49 


71.14 


0.0 


420 


10 QUANTILE 


79.37 


79.04 


78.33 


74 . 


25 


82.87 


82.52 


81.89 


88 . 


50 


S6.29 


85.89 


85.42 


100. 


75 


89.50 


90.12 


88.69 


113 


90 


92 . 50 


92.07 


91.56 


124 


I .XI MUM 


104 .16 


103.61 


112.79 


163 



4 

6 

00 

1 

20 

24 

17 

2 

5 

0 
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CASE VI. WEIBULL DISTRIBUTION w = 1.5 

2 

Y = (- 1442695.04 In U) 3 



MOE 


CEP 1 


CEP 2 


CEP 3 


CEP 4 


MEAN 


89.58 


89.28 


87.51 


102 


STD DEV 


7.07 


7.04 


12.85 


25.7 


MSE 


158.65 


164.46 


321.29 


664 


GED MEAN 


89.30 


89.00 




98.4 


VARIANCE 


50.07 


49.54 


165.09 


660 


COEF VAR 


0.08 


0.08 


0.148 


- 0.25 


SKEWNESS 


0.07 


0.07 


- 4.27 


- 0.04 


KURTOSIS 


0.22 


0.21 


27.91 


0.11 


RANGE 


44 . 55 


44.26 


125.57 


161 


MINIMUM 


69.36 


69.10 


0.0 


31.5 


.10 QUANTILE 


80.27 


79.99 


78.95 


67.6 


.25 


85.39 


85.12 


83.80 


86.4 


o 

in 

• 


89.59 


89 . 25 


88.60 


101 


.75 


94.37 


94.07 


93.44 


120 


.90 


98.31 


98.00 


97.48 


133 


MAXIMUM 


113.91 


113.37 


125.57 


19 3 
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CASE VII 



WE I BULL DISTRIBUTION w= .5 



MOE 


CEP 1 


MEAN 


166.57 


STD DEV 


42.61 


MSE 


6246.94 


GED MEAN 


161.26 


VARIANCE 


1815.38 


COEF VAR 


0.26 


SKEWNESS 


0.62 


KURTOSIS 


0.53 


RANGE 


250.79 


MINIMUM 


75.81 


.10 QUANTILE 


116.20 


.25 


135.54 


.50 


163.59 


.75 


191.58 


V£> 

O 


220.85 


MAXIMUM 


326.60 



In U) 2 



A A 



2 


CEP 3 


CEP 4 


05 


87.45 


146 


99 


31.55 


99.2 


96 


1153.15 


11956 


32 


81.72 


113 


96 


995.65 


9840 


26 


0.36 


0.68 


70 


0.76 


1.39 


76 


1.48 


3.26 


53 


218.72 


746 


48 


18.67 


3.24 


74 


49.30 


41.2 


27 


66.56 


75.8 


51 


84.41 


126 


45 


105.80 


194 


60 


127.69 


282 


01 


237.40 


750 



4.27 

CEP 

170 

44 

6930 

164 

2023 

0 

0 

0 

258 

78 

117 

137 

165 

196 

227 

337 
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CASE VIII. WE I BULL DISTRIBUTION w= .1 
Y = (- 3.623886 In U) 10 



MOE 


CEP 


1 


CEP 2 


CEP 3 


CEP 4 


MEAN 


2.6 


10 5 


3.6 10 5 


131.06 


5.1 10 4 


STD DEV 


7.8 


10 5 


1.1 10 6 


335.02 


2.8 10 5 


MSE 


6.7 


10 11 


1.4 10 12 


1.13 10 5 


8.2 10 10 


GED MEAN 


6.0 


io 4 


7.4 10 4 


85.89 


674 


VARIANCE 


6.0 


10 11 


12 

1.3 10 ^ 


1.1 io 5 


7.9 IO 10 


COEF VAR 


2.97 


3.14 


2.56 


5.43 


SKEWNESS 


7.95 




8.38 


14.16 


14.4 


KURTOSIS 


81.28 


89.17 


232.51 


255 


RANGE 


9.8 


10 


1.4 10 7 


6214.98 


5.4 10 6 


MINIMUM 


461. 


12 


476.59 


5.65 


7.2 10 6 


.10 QUANTILE 


6859 


.66 


7661.19 


30.28 


4.41 


in 

CM 

• 


1.8 


10 4 


2.2 10 4 


54.43 


77.0 


o 

in 

« 


5.9 


10 4 


7.4 10 4 


93.37 


1070 


.75 


1.7 


10 5 


2.2 10 5 


140.28 


9280 


.90 


5.6 


io 5 


7.3 10 5 


195.25 


8.3 10 4 


MAXIMUM 


9.8 


io 6 


1.4 10 7 


6220.63 


5.4 10 6 
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COMPARISON OF ESTIMATORS 



USING MEAN SQUARED ERROR AS A CRITERION: 
MSE [CEP ] = O ~ + [E ( CEP) - CEP ] 2 



-N A- V 



CASE NO. 


DISTRIBUTION 




CEP 1 


CEP 2 


CEP 3 


CEP 4 


I 


RAYLEIGH 






136 


136 


165 


1549 


II 


UNIFORM 






109 


107 


368 


1744 


III 


POWER 






2460 


83,000 


444 


2279 


IV 


TRANSLATED 


RAY. 


134 


134 


177 


1524 


V 


WEIBULL w 


= 


2. 


218 


228 


379 


382 


VI 


WE I BULL W 


= 


1.5 


158 


164 


321 


664 


VII 


WEIBULL w 


=: 


.5 


6246 


6930 


1153 


11956 


VIII 


WEIBULL w 




.1 


6.77 10 11 


6.78 10 11 


1.13 1C? 


8.17] 
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VII. CONCLUSION 



In general, we found the maximum likelihood power 
estimator cfiP 3 is better than moment method power estimator, 
in that it does not have the limitations and exists for all 
parent distributions we considered. Using minimum mean 
squared error as a criterion, the estimators based on the 
power model have inferior performance in the "shallow tailed" 
cases, but were best in the heavy tailed cases, this is 
especially true for the w = .1 Weibull case, in which it 
has less bias than its competitors and has relatively small 
variance . 

We found the Rayleigh based CEP estimators (CEP 1) to 
be possibly the best overall, especially for the shallow 
tailed cases. The maximum likelihood jackknife Rayleigh 

A 

estimator (CEP 2) had no better performance than the unbiased 
Rayleigh estimator. 

A 

The value of the sample median estimator (CEP 4) can 
be found for small samples with such rapidity that in certain 
cases the time saved may compensate for the accuracy lost. 

In practice, when analysts find several impact points 
comparatively far away from target, they should choose 

A 

maximum likelihood power estimator (CEP 3) . Otherwise they 

A 

could choose Rayleigh based CEP estimator (CEP 1) . 
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APPENDIX 1 



MAXIMUM LIKELIHOOD POWER ESTIMATOR 



DIMENSION V ( ZOO) , S (2001), A (2001) 

DIMENSION CEPI5G0) 

CALL OVFLCW 

I S EED= 4 70 7 23 

READ (5,100) AO, SO 

READ (5,200 NS, IST0P,S(1) , A ( 1) 

S N = F L 3 A 7 (NS) 

F A =— 1 » 0 /A 0 
KS=500 

DC 1 K I — 1 , KS 
C 

C GENERATE power SQUARE deviate 
DC 5 i\ = 1 , N S 

V ( N ) = C . C 

CALL RANDOM ( IS H ED. U, 1 ) 

Y ( N ) = SO *( U* * F A - 1 .0 I 
5 CONTINUE 

C 

C SOLVE NONLINEAR EQUATION UNDER NEWTON'S SECOND METHOD 

DO 10 1 = 1, I STOP 
I T tRS = I 
A ( I + 1 ) - 0. 0 
S( I+l)=C.O 

CALL SF (SC 1 ,Y,NS,B) 

CALL CF ( S ( T ) , Y , H S , C ) 

CALL DF (SC ) ,Y,NS,0) 

E = A L 0 G ( S< I ) ) 

UX=5N/A ( I )+SN*E-E 

VX=SN* A (I )/S (I ) — ( A ( I ) + 1.0)4D 

WX=-SN/(A( I ) *v2 ) 

X X = S N / S (D-D 

ZX = -S N ;; A ( I )/ (S ( I )**2 ) +( M I ) +1.0) ‘AC 
DET=WX*ZX-XX**2 
C 

C CHECK THE OE T F 0 M 1 N A T I ON 

IF (C5T.EQ.5.0) CO TO a? 

UU=ZX*UX-XX~VX 

VV=WX*VX-XX*UX 

C 

C WHETHER TO STOP "HE ITERATION 

C H T = - ( LX*U'J-rVX*VV ) /DET 
IF (CHI .LE. ) • C 4 0 ) GO TO 03 

A ( 1 + 1 ) = A ( I J-'JU/YFT 

S ( 1 + 1 ) = S( I )-VV/DcT 
IF (S(I+1 ) .GT.G.Q) GO TO 10 
S ( I + 1 ) = S ( I 1/2.0 

10 continue 

3 £ AOUT = A( ITERS ) 

SGUT= S ( ITERS ) 

C 

C ESTIMATE CEP 

AK = 2**( 1. O/AOUT ) -1.0 

CEP ( KI ) = (S3UT*AK)*=?G.5 

WRITE ( 6 , 40 J ) SOUT, AOUT, ITERS, CEP(KI) 

1 CONTINUE 

CALL HISTG (CEP,KS,0) 

99 STOP 
END 



SUE ROUT IN E FF <S,Y.NS,B) 
DIMENSION Y ( 2GG ) 

B = C.O 

DC 30 N= 1 , NS 
B=B+ALOG(S+Y(N) ) 

30 CONTINUE 
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RETURN 

END 

SUBROUTINE CF (S,Y,NS,C> 
dimension y(200) 
c = o.o 

DO 40 N - 1 * NS 
C=C+1. 3/(S+Y(N) ) 2 
40 CONTINUE 
RETURN 
ENC 

SUBROUTINE OF <S»Y,N$,D> 
DZVENSICN Y( 200 ) 

0 = 0.0 

DO 50 N=1,NS 
D=D+I.O/(S+Y(N) ) 

50 CONTINUE 
RETURN 
END 
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APPENDIX 2: RAYLEIGH UNP-IASED AND JACKKNIFE ESTIMATOR 



D I MENS ION Y{ 50, 50) ,SP< 50 J , PSE ( 50 ) , CEP( 500) 

DIMENSION RAY (500) 

C 

C GENERATE RAYLEIGH DEVIATE 

CALL OVFLOW 
I ScED=47C723 
C 

C READ IN ORIGINAL POP. PARAMETER VC , TOTAL SAMPLE SIZE 

C NS, SU6GRCUP f.AMPL c SIZE KS 

READ (5,100) VO , NS , KS 
CPC*- 1. 17741G02*SQRT( VO ) 

MS=N5/KS 
SM=FLOAT(MS ) 

SN=FL3AT ( NS ) 

SK=FLGA7 ( KS ) 

WRITE ( o , 20J ) VO, CPC, NS, MS , KS 
K E = 5C 0 

DO 1 KR E = 1 , KS 
DO 5 M = 1 , MS 
DO 5 K= 1 , KS 
Y ( M , K ) = 0 . 0 

CALL RANDOM (ISEED, U, 1 ) 

Y(M,K)=1C000.0*( l.D/U-1.0) 

5 CONTINUE 
C 

C ESTIMATE THE PARAMETER BASED ON ALL DATA 

T = C.O 

DO 10 M = 1 , MS 
CO 10 K = 1 , KS. 

T = T +Y ( M , K ) 

1C CONTINUE 

T M = S C R T ( T / ( S N - 2 • 0 ) ) 

C 

C MAXIMUM LIKELIHOOD UNBIASED ESTIMATES 

R AY ( KRE )= 1.1848*TM 
C 

C ESTIMATE THE PARAMETER BASED ON ALL DATA EXCEPT THOSE 

C IN THE ITH GPGU? 

DO 20 M= 1 , MS 
S = 0.0 

DC 25 1=1, MS 
IF (l.EQ.M) GO TO 25 
DC 26 K = 1 , K S 
S = S + Y ( I ,K) 

26 CCNTINJE 
25 CONTINUE 

SP (M ) - SQKT ( S/ ( (SN-SK) *2.0) ) 

20 CONTINUE 
C 

C COMPUTE THE MS PSEUDO VALJES 

DO AC m=i,ms 

P S E ( M ) = 3 M * 7 M - ( 3 M - 1 . 0 ) * S P ( M ) 

AO CONTINUE 
C 

C ESTIMATE JACKKNIFE PARAMETER 

ACK=C .0 
DO 45 M=1 , MS 
AC K = ACK + PS E ( M) 

45 CONTINUE 

VJACK=ACK/SM 
S C K = 0 .0 

c 

C COMPUTER THE CONFIDENCE INTERVAL 

00 50 1=1, MS 

SCK=SCK+( PSE ( I )-VJACK)**2 
5G CONTINUE 
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SVJA=SCK/ ( .0) ) 

U P E = V J A CK-2. 0*SQ.RT (SVJM 
PL B =V J AC K +2 • O^SQKT (SVJ4) 

C 

C ESTIMATE C5? FROM JACKKNIFE PARAMETER 

C E P ( K R F l=1.1774*VJACK 
WRITE (6,303) V J AC K , UP B , PL 6 , CE P ( Kk E ) 

1 CONTINUE 

CALL HISTG <RAY,KE,C) 

CALL HISTG (CEP,Kf,G) 

STOP 

END 
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